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Abstract 

We derive a model Hamiltonian whose ground state expectation 
value of any two-body operator coincides with that obtained with the 
Jastrow correlated wave function of the many-body Fermi system. 
Using this Hamiltonian we show that the variational principle can be 
extended to treat systems with dynamical mesons, even if in this case 
the concept of wave function looses its meaning. 
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Variational methods, based on the Fermi Hyper Netted Chain (FHNC) 
expansion |jl], ||, |^ , also connected with the Correlated Basis Function pertur- 
bation theory 0, have in the last years provided a satisfactory description 
of the static properties and of the response functions of the nuclear matter 
and of other quantum fluids like, e.g., the liquid ^He. 

The FHNC approach is however seemingly constrained to non relativistic 
potential theories, where the wave function concept makes sense and corre- 
lated wave functions of the type: 

= G'0(n,...,rA), (1) 

with 0(ri, . . . , ta) being a Slater determinant and G being a correlation op- 
erator, can be used. The simplest of such correlation operators is given by 
the Jastrow Ansatz ^ 

Gi = n/M- (2) 

In this letter we show that a variational method exists that is able to properly 
account for systems where the wave function concept breaks down: this 
happens for instance when dealing with dynamical meson exchange, a process 
not reducible to the frame of a potential theory, which will be the main topics 
of the present letter. The same underlying idea however could be applied to 
other fields in solid state physics (Ising and Hubbard models, for instance), 
and could provide further improvements in the customary FHNC calculations 
since our approach, as we shall see, can allow an energy-dependence in the 
correlation function. 

In this letter we wish to outline the general idea, namely the translation 
of the variational problem from the setting of a trial wave function to the 
setting of a trial Hamiltonian. Thus we are able to overcome the constraint 
represented by the need of defining a wave function, so allowing the treatment 
of dynamical meson exchanges (or similar situations). We exemplify this 
treatment for the case of the Jastrow correlated wave function. 

We avoid here technical difficulties by confining ourselves to nonrelativis- 
tic nucleons coupled to mesons but disregarding for instance state-dependent 
correlations; more detailed problem will be dealt with separately in subse- 
quent papers. 

To achieve our goal, we have to solve three problems: 
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1. to embody the variational principle into a formalism able to handle 
effective Hamiltonian, [0, || 

2. to translate the variational wave function (we will limit ourselvs to the 
Jastrow Ansatz) into a Hamiltonian formalism, 

3. to find an effective Hamiltonian with nucleons only as true degrees of 



freedom but exactly including the dynamical meson exchanged, 10, 11 



These items are nicely expressed in the Feynman Path Integral language 
[0]. Concerning item 1 the variational procedure is derived as follows. Let 
H[ip\ Tp) be the Hamiltonian of a many-body system and 

-/ dTH{ip^,ip) 

V[i,\i,]e , (3) 

the partition function. The ground energy of the system is 

Bo = -l..ni^. (4) 

The partition function is of course defined in the euclidean world, so that the 
Feynman-Kag measure and consequently the integral in (j^) are well defined. 
We introduce a weighted average on the space of the Feynman-Kag- 



integrable functions as 



- / drHo 

< J >Ho= ^ , (5) 

- / drHo 

p 

- J drHo 

where e » plays the role of a weight function. Actually Hq will be 
interpreted as a suitably chosen model Hamiltonian either solvable or at 
least easily manageable. It is worth stressing that Hq is not necessarily 
the representation in the Bargmann-Fock space of a true Hamiltonian (i.e. 
an energy-independent hermitian operator bounded from below): it is only 
requested to be real and positive-definite in the euclidean world. Thus an 
explicit frequency dependence of Hq can be allowed. 
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Using d^) and the inequality < e ^ >> e '^^^ we find 

-J{H-Ho)dT . -JdrHo 

Z = <e j V[^\ij]e 


-< J iH-Ho)dT>HQ /. - / drHo 

> e / V[tp^,^]e , (6) 

from wliicli tlie result 

< JdrH >Ho 

Eo<E,ar = }ira , (7) 

easily follows. 

As for Ho, H is not requested to be hermitian or energy- independent: it 
only needs to be real in the euclidean world. Thus (0) applies to effective 
Hamiltonians as well, provided they are real after a Wick rotation. 

This procedure can be applied to a meson theory provided an effective 
Hamiltonian is derived; on the other hand we are free to choose a suitable 
Hq. This makes the difference with the variational approach: there we guess 
the ground state wave function, here we guess an unperturbed Hamiltonian. 
How could we force the two approaches to coincide? 

Now we are ready to better formulate question 2): we look for an operator 
Ho such that the Jastrow wave function F will be its ground state. Presently 
we need something less, namely that E^ar coincides with the expectation 
value of H on the Jastrow wave function ipj of eqs. 

Let us consider the following model Hamiltonian: 

^o"'^ = ^J2 + Y ^ K(lWk+galap_qap . (8) 

k<kp kpq 



We claim that 

_ < ^j\H\i,j > 

var ^ I \ I ^ ' 

provided the correct link between A(g) and f{r 12) is found. 



(9) 
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We will present here only the proof that the above equality holds for the 
expectation value of the potential energy, namely 



/3 

< / drV >Ho / ; IT/I / ^ 

e> = ii„, . . (10) 

Since this result will be true independently from the particular V present 
in the Hamiltonian, it implies that the diagonal part of the two-body den- 
sity matrix either derived from our model Hamiltonian or from the Jastrow 
Ansatz will coincide. Consequently the one-body density matrix will do the 
same and the kinetic energies will also coincide. 

We shall show that Vj^J is given by the sum of the FHNC diagrams. 

Let us expand the Ihs of eq. (|^) in Goldstone diagrams using the bilinear 
part of Hq"^^ as unperturbed Hamiltonian. A simple analysis shows that at 
each order the quantity V^^J is independent from A (a n^^ order diagram has 
n AA(g) terms and n energy denominators: since each of the latter oc A, 
then all A's cancel). On the contrary the Green's functions of the theory 
depends upon A and becomes singular when A — >■ 0. 

Let us now show the connection between the perturbative expansion of 
Vj^J and the FHNC diagrammatical expansion of < '?/'j|V^|'?/'j > / < ipj\ipj >. 

We consider a Feynman diagram of the expansion of < V >Ho {s^ first- 
order diagram is shown in fig. la). The solid lines (unperturbed Green's 
functions) are given by 

Ko + irj Kq — A — IT] 

We split Qo into two pieces: 

go{k,ko) = go{k,ko)+K{k,ko) (12) 

^°^^'^°) - '^'^-'^{k;^x^,-k^} (1^) 

Next we sort new diagrams from the Feynman graph by expanding it by 
means of (|12|). In our example three graphs are generated (fig. lb): single 
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lines now represent go lines and double lines K lines. Next we carry out 
the frequency integration, that will be in general elementary if A is energy- 
independent, so coming to Goldstone diagrams. In each of them the K lines 
(which are momentum-independent) cannot leave any trace and are shrinked 
to points. Tadpoles also translate into constants and are also represented as 
points. The result is shown in fig. Ic. 

All the FHNC diagrams can be reproduced in a similar way. In fact the 
frequency integration reduces the product of the energy denominators in a 
given diagram to a factor niA^ (with m integer, k being the order of the 
diagram), and this factor is canceled by the coming from the interaction 
terms. Thus the K line can only influence the integer m. The go lines after 
the frequency integration become ^-functions. We remind the reader that 
the Pauli correlation line of the FHNC diagrams corresponds to the Fourier 
transform of the ^-function, i.e. 

i{kpx) = -^ J ^^e{k, - pY"^- . (15) 

The multiplicity factors v coming from spin-isospin summations necessarily 
coincide both in the present approach and in the FHNC formalism, since they 
come from the counting of closed loops; the powers of the density, coming 
from integrations over single go lines, also coincide with those of FHNC for 
dimensional reasons. 

So far we have demonstrated the topological equivalence between the 
FHNC diagrams and those of the present theory. It remains clarify the 
relation between the correlation function h = p — 1 oi FHNC theory and the 
"potential" A and then to verify that the coefficients carried by the diagrams 
coincide in the two expansions. 

For these purposes we follow the derivation of FHNC diagrams given in 
. There the two-body operator U = 2 log p is introduced and a fictitious 



time dependence is assigned both to U and to the potential V with the results 

< $o|1/e^|$o >=< %\T /V(t)e/o ^W-^^l^o > . (16) 

Jo 



Next, in ref.|13| the r.h.s. of the above equation has been expanded by means 
of the Wick theorem, and the FHNC diagrams originally derived in ref. 
are recovered. 
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Now, within the present approach, let us consider the Feynman-Kag in- 
tegral 



/ 



/5 

■ J H^/\r)dT .(3 



P[a,(r),4(r)]e o / V{Od^ , (17) 

Jo 



for our infinite system: since the "unperturbed part" of the model Hamil- 
tonian has the Fermi gas Slater determinant |$o > as a ground state as in 
(|16D, we can reinterpret (|l^) as 

<%\T r V{Od^e-^fo^^^^'^^\%> . (18) 
Jo 

Since the physical quantities we are interested in are independent of A, we 
can put A = and make the replacement r — r//3 to recover (|1^) provided 
we identify U with —A in configuration space. Next we take the limit /3 — > cx) 
(and consequently A — ^ 0) in order to apply (H) , (^ and (^. Since the value 
of A is immaterial, the equivalence between the FHNC development and the 
perturbative expansion of our model Hamiltonian is proved. A similar proof 
can be done for the kinetic energy part of the Hamiltonian. 



Finally, in |T3|, it is shown that the ladder series of fig. 2 can be summed 
explicitly at the prize of redefining U. The diagram structure being identical, 
we can do the same in our case by means of the identification 

h = e-^ -1 , (19) 

which, we remind, is justified in the limit A — > 0. 

This completes the second step of our program. Finally, to apply the 
previous results to pion physicis, we need an effective Hamiltonian containing 
fermions only, but also accounting for pion exchanges. In the context of QFT 
an effective action has been derived in ref . 0] . Following the same path we 
can write a lagrangian density of the form 

C = Cn + + i'^\x)r{x)'^{x)(f){x) (20) 

Cn and £^ being the free lagrangians for nucleon and pion fields respectively 
and r a spin-isospin operator (plus eventually derivatives) and then the cor- 
responding generating functional (the partition function in the Minkowski 
world), here for simplicity with fermionic sources only: 

Z[r],r]^] = /p[^,7/'^0]e'/'^"[^+'^'^+^^'^] (21) 
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Since the pionic field is at most quadratic in the exponent, the integral over 
is gaussian and can be carried out explicitly, getting 

where 

-S'eff = J dxdy^\x)GQ^{x - y)'ip{y) (23) 

~\ I (^^dyilj\x)r{x)ilj{x)Do{x - y)i!^{y)r{y)i!{y), 

Go being the free fermion propagator and Dq the (fully dynamic) pion prop- 
agator 

'''' 

The last step is to apply eq. (|) and (0) to the action (^). Thus first 
we go to the euclidean world, so that Dq becomes real and bounded from 
below. Next we can use (^, since the qo dependence in Dq does not prevent 
eq. (|^) from holding, and we expand diagrammatically the average value of 
the (euclidean) action Ses using Hq^^ as model Hamiltonian. 

The diagrams so obtained have the same topological structure as the 
standard FHNC diagrams previously discussed, provided the time variable 
is added for the external points. Thus these diagrams can also be summed 
up by means of the usual FHNC technology, but now the constant A is no 
longer irrelevant: in fact in one energy denominator (namely the one of the 
meson) it is compared with a finite quantity (the meson mass). Thus A can 
be regarded as an extra variational parameter. 

A more realistic variational Hamiltonian than Hq^^ should be used in 
connection with effective Hamiltonians of the type discussed above. The 
interaction term in Hq"^^ should be modified to include state-dependent cor- 
relation effects present in a realistic Moreover, a time (or frequency) 
dependence could be introduced in A, which would imply adding a time de- 
pendence in the internal points of the FHNC diagrams, or in other words to 
consider time dependent correlations. Work in this direction is in progress. 
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Figure Caption 

1. la A first order (with respect to A) Feynman diagram, 
lb The diagrams obtained from la by means of eq. (|T2p. 
Ic The corresponding FHNC diagrams. 

2. The ladder series corresponding to eq. (|l^) 
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